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Abstract 

We extend the problem of (1 + 1) circular dilaton gravity to include charged 
particles. We examine the two (charged) particle case in detail and find an 
exact equilibrium solution. We then extend this to N— particles and obtain a 
solution for this case as well. This class of solutions corresponds to TV-particles of 
the same mass, spaced evenly around the circle with charges chosen so that the 
electric field satisfies E 2 =constant. We discuss the relation of these solutions to 
the previous uncharged equilibrium solutions and examine the behavior when 
the number of particles is large. We comment on the challenges in further 
generalizing the solutions we obtain. 

1 Introduction 

Circular gravity is a general relativistic gravitational system in (1 + 1) dimen- 
sions in which the topology of space is equivalent to that of a circle S . Its 
counterpart, lineal gravity (the same system with the S 1 topology replaced by 
M 1 ), has been the subject of much investigation over the past twenty years, mo- 
tivated by the desire to understand black holes pj], space-time structure pUE], 
canonical quantum gravity string physics [Sj, information loss in black hole 
evaporation |Zj and the N— body problem [HIE]. 

This latter problem involves consideration of a collection of N particles 
interacting through their own mutual gravitational attraction, along with other 
specified forces. Such systems are referred to as self-gravitating systems. Apart 
from functioning as prototypes for higher dimensional gravity, in lineal gravity 
they also approximate the behavior of some physical systems in three spatial 
dimensions, such as stellar dynamics orthogonal to a galactic plane, the collisions 
of flat parallel domain walls, and the dynamics of cosmic strings. The ergodic 
and equipartition properties of non-relativistic one-dimensional self-gravitating 
systems are a current subject of research an d recently such systems were 
shown to exhibit fractal behavior 

In the past few years these systems have been extended to relativistic lineal 
gravity. The 2-body problem has been solved exactly ^2], as have its extensions 
to include cosmological expansion ^J^] and/or electromagnetic interactions 
|15j . These extensions, while considerably more complicated, are nevertheless 
amenable to statistical-mechanical analysis (§] . They include the non-relativistic 
self-gravitating systems ^01 D3| as a limiting case. 

Comparatively little is known about the circular counterparts of these lineal 
systems. Preliminary investigations |17| have shown that qualitatively new fea- 
tures emerge that are absent in the non-compact lineal case. For comparison, 
in the absence of other fields there are no static non-relativistic equilibrium so- 
lutions for N particles on a circle. This is because the non-relativistic potential 
grows linearly with increasing distance from the source, and a non-relativistic 
circular topology admits no solutions to the matching conditions for any N > 1. 
However it has recently been shown that the relativistic case has an exact equi- 
librium solution for N particles £Q . This solution is possible because the space- 
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time itself either expands or contracts in response to the sources, an option not 
available in the non-relativistic case. 

In this paper we extend the formalism of the iV-body problem in circular 
gravity to include electromagnetic interactions. We search for equilibrium so- 
lutions in the case when the particles are charged. We examine the charged 
2-body system in detail, and find a class of charged solutions that can be gener- 
alized to ./V-particles and encompasses the previous non-charged solutions. Our 
solutions are not the most general, and we comment on the challenges involved 
in extending the solutions that we obtain. 



2 Canonical Circular Gravity 

We begin with the general framework previously formulated for canonical cir- 
cular gravity minimally coupled to N particles in the presence of a cosmological 
constant A [TJ\, extending the formalism to include coupling to electromag- 
nctism. This can be done in a manner analogous to that way that charges were 
added to lineal gravity JS]- The resulting action is 
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where ^> is the dilaton field, and F^ v are the vector potential and the field 
strength, g^ and g are the metric and its determinant, R is the Ricci scalar, 
and T a ,z a , and e a are the proper time, position, and charge of the a-th particle, 
with k — 8irG/c 4 . The range of the spatial coordinate x is — L < x < L, and 
the circular topology implies that all fields must be at least C 1 functions of x 
with period 2L, therefore 



f(L) = f(-L) and f'(L) = f'(-L) . (2) 

The action describes a generally covariant self-gravitating system (with- 
out collisional terms, so that the bodies pass through each other); its field 
equations are 
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(3) 

kT^ + -g^A , (4) 

(5) 

(6) 
(7) 

where the stress-energy due to the point masses and the electric field is 

Tp, v = X! TOa [ dT a^=9^9u P ^ ± ^f^S 2 (x - Z a {T a )) (8) 

a J V 9 dT a (lT a 

+ — g {F^F v ^ - i g ^F a pF al3 } : 

Recall that there is no magnetic component of the field in (1 + 1) dimensions. 
Conservation of T M „ is ensured by eq. (@J . Inserting the trace of eq. Q into 
eq. J2J yields 
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d f d<\ i .d^dz*" 



R-A = K T^. (9) 

Consequently eqs. ©-0 and © form a closed system of equations for 
gravity matter, and electromagnetism leaving (@J) to determine the evolution 
of the dilaton field '5. The preceding system generalizes Jackiw-Teitelboim 
lineal gravity |18) . which equates scalar curvature to a (cosmological) constant 
A, reducing to it if the stress-energy J3J vanishes. The specific choice of dilaton 
coupling in R = T theory is such that the evolution of the dilaton does not 
modify the aforementioned reciprocal gravity/matter/electromagnetic dynam- 
ics. The total system is a closed system that describes the evolution of the 
single metric degree of freedom and the N degrees of freedom of the charged 
point masses. The evolution of the charged point-masses governs the evolution 
of the dilaton via 10} , whose divergencelessness is consistent with the conser- 
vation of T^ v) yielding only one independent equation to determine the single 
degree of freedom of the dilaton. 

The metric is written as ds 2 = — {N^dt) +7 ( dx + ^dt) , where No is the 
lapse and N\ is the shift. With this metric the extrinsic curvature K is then 

K = (2N j)- 1 (2d 1 N 1 - 7 - 1 iV 1 a l7 - d j) = ^(tt - II/7) (10) 
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where the conjugate momenta to 7 and ^ are respectively denoted ir and II. 
We can take the extrinsic curvature K to be a time coordinate r (t) of the circle, 
thereby allowing the elimination of ir from all field equations since ir = + ~ . 

Simplifying the action JTJ involves decomposing the scalar curvature in terms 
of the extrinsic curvature via the equation ^/—gR = —2do(y/jK) + 2di[{N\K — 
diNo)/y/j. Substituting this result into the action and eliminating 7r via 

7T = ^ + S yields [El H3 

d 2 x^2p a z a S(x^z a )+n^- t ^ + ln 7 )-n^ , (11) 

where H = <f dxTi. =— § dx^j. H is the circumference functional of the circle 
when f is constant. The dynamics is that of a time-dependent system, with the 
time-dependence corresponding to the time-varying circumference of the circle 
of constant mean extrinsic curvature. As with general relativity on spatially 
compact manifolds in (2 + 1) dimensions |19j . the spacetime metric induces the 
metric da 2 = gn (x, t) dx 2 on each t ^constant hypersurface. Since any non- 
singular metric on a circle is globally conformal to a flat metric, we choose the 
spatial metric to be such that 7 = 7(f) , i.e. 7' = 0. In other words, we write 
the gn component of the metric as "/(t)h(x,t), and then set h(x,t) = 1 by a 
coordinate transformation. By choosing a time parametrization so that f^/7 is 
constant, the Hamiltonian can be made time-independent. 

It is important to note that one of the intermediate steps in the previous 
simplification is a constraint equation on the electric field E |15|. which is 

^- = J2^ S (x-z a (t)) (12) 

a 

The solution to this constraint equation (|12|) is 

E=-J2e a e(x-z a (t))+E (13) 

a 

where e(x) is the step function e(x) — and e(0) = 0, or more explicitly 

e(a;) = —l,x < 
e(x) = 0,x = 
e(x) = l,x > 

We chose a regularization such that f(z a ) = l/2(f(z a + e) + f(z a — e)) for 
any quantity that has a jump discontinuity. This is necessary as otherwise 
it would not be possible to satisfy field equations that depend on the specific 
value of the fields at the location of the particles (i.e. field equations Q2() [l .lj21 [l 
or cquivalently (1291 1. (|3()ll . ). The choice of e = at the particle is a reflection 
of this regularization and is necessary to form a solution. 



4 



Applying the smoothness condition across the boundary (J2J to the equation 
for the electric field (|13|) implies that J2 a e a = 0- So the total charge of the 
(1 + 1) circular gravity system must equal zero. 

The resulting field equations from the action (|llfl therefore become 

2W -IW' + Y,Pa6(x- z a {x )) = (14) 

a 

*" - \^') 2 - W 2 + 7 ( r2 - A / 2 + ^ 2 /2) 

+ «EVpI+^<5(x-z a ( 2 ; )) = 0, (15) 

a 

Aft - f 7 

-aJt^-A^ + kSV^ + kV iml 6(x-z a (x°)))=0, (16) 
A{ - 7 /2 + 7 rA = 0, (17) 

rt + 9 l( -iiv 1 n +5 ^_W + -l f ^) = o, (is) 

<i- + 2A (k-^+t) -Ni (-#') = 0, (19) 
V V7 / 7 



. + cWo /pg | m2 N pj c>7 <9Ai p a | ^ p a d 7 

dz a y 7 a 2^^. + TO 2 7 2 <9z a 7 1 J 2 dz a 

dxN ^yE— = 0, (20) 

A r=?= + — = . (21) 



In the equations l|20|l . (|21|l the prime denotes derivative with respect to x 
and the dot denotes derivative with respect to t. All metric components (Nq, 
Ni, 7) are evaluated at the point x = z a where 



df _ df(x) 



dz a dx 



(22) 



Previous work in the e Q = case £Q has shown that an equilibrium solution 
exists when all particles have the same mass and are evenly spaced around the 
circle. So it is reasonable to search for (possibly similar) equilibrium solutions 
in the charged case. Introducing the ansatz m a — m,z a = 0,p a = into eqs. 
<|14l) - (|21[) leads to the simplification 
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211' - II*' = 0; (23) 

*" - j(* ) 2 - («n) 2 + 7 (j 2 g — 

+K^2y^mS{x - z a ) = 0(24) 

Aff - T7 - JV j 7 f r 2 - A ~ )+k/2\] ^mS(x -z a )\ = 0(25) 



a ) 



N[-j/2 + -/TN = (26) 
n + d^U + ^N^ + ^N') = 0(27) 

i> + 2N Q ( K^L+t| -JVi(i*') = 0(28) 
V V7 / 7 

dN f , Ar 



d ^-m + J dx N ^E- = 0(29) 

= 0.(30) 

3 Two Neutral Particles of Unequal Mass 

We begin by attempting to solve the case of two uncharged bodies of unequal 
mass. Here we find that the resulting constraint equations seem to force very 
specific conditions on the system. We shall encounter a similar situation in the 
charged case 

As before we use the same field equations l[14 jl -i|21 [l and again equilibrium 
solutions are characterized by z a — = p a , corresponding to a situation in which 
all particles arc motionless at various points around the circle. This implies from 

mrm that 

f^=0 N 1 (z a ) = . (31) 

the field equations p4H19(l can be solved in terms of the quantity c 2 = 7 (t 2 — A/2) , 
which is a function only of the time coordinate. 
The solution to ( 114ft is simply 

n = n (i)exp*/2 (32) 
which can then be substituted into 1|15|) to obtain 

JV 

[ty') 2 - (VrTU 2 fixn * 4- r 2 + 



1 

*" " t(*') 2 " («n ) 2 exp * + c 2 + K Vl m a S(x - z a ) = (33) 



a=l 

for a system of N particles in equilibrium on the circle. 
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When N = 2, the two particles are placed on a line at coordinates (z\, z 2 ) 
respectively, dividing the line into three regions. Without loss of generality, a 
shift in the x-coordinate allows us to place the origin equidistant from either 
mass, so that z\ = z = —z 2 . 

In the 2-body case the matching conditions imply from (|33[l that 

= * 3 (-L) & 1 (L) = %(-L) (34) 

= * 2 («) *s(-*) = *2(-«) (35) 

^ 2 (z) = -/«V7mi « / 2 (-*)-*3(-«) = -Ky/jm 2 (36) 

for the dilaton field. The conditions l|34|l are smoothness conditions at the join 
of the outer regions where the line segment is identified to a circle. Conditions 
(|35(l are continuity conditions at the locations of the particles, and conditions 
l|36l) arc junction conditions which follow upon integration of l|33|) . 

In general there are three possible choices for c 2 : c 2 = 0, c 2 = cP\_ > 0, 
c 2 = cL < 0. It is straightforward to show that no solutions exist for c 2 = 0. 
The solution in each region 1=1,.., 3 when c 2 = is 

i = -2 In [«n (Bj + ( Sl x + (1-2) L))] (37) 

where B\ are functions only of time. The constant (1-2) L is inserted for con- 
venience, and s 2 = 1 for 1=2 and sj = 1 otherwise, eq. (I34|l implies that 
B\ = Bs = B, and (|35fl is easily shown to imply that B 2 = B, and z = L/2, 
with S2 = — 1; the only alternative solution yields L = 0, which is unacceptable. 
Finally, equation (|36|l implies after some algebra that 

1 1\ . -4 / I J_ 



B =_ L=^ — + — (38) 

which is not allowed, since L > 0. Hence there are no equilibrium solutions if 
c 2 = 0, consistent with the single-particle case. 

Consequently we need only examine the cases c 2 = c 2 ^ > and c 2 = c?_ < 0. 
Our analysis shall concentrate primarily on the former case, since the latter case 
can be obtained from it by setting c — ► ic in our solutions, with only slight 
modifications that we shall discuss as they arise. 

For the case c 2 = c 2 , > 0, the solution in each region 1=1,.., 3 is 



#i = -2 In [Bj - cosh (c+ (x + (1-2) L) + Aj)] - 2 In 



(39) 



where as before (A\,B\) are functions only of time, and the constant (1-2) L 
is inserted for convenience. Smoothness at the identification point yields from 
(3D that 

Si - cosh (Ai) _ B 3 - cosh (A 3 ) 



sinh(Ai) sinh(A 3 ) 



Bj - cosh (Ai) B 3 - cosh (A3) 



(41) 
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It is straightforward to show that the only solution to these equations is A\ = 
A 3 = A,B 1 = B 3 = B. 

Consider next matching the solutions at the points (21,^2)- Eqs. i|35|) yield 



B - cosh (c + (z - L) + A) 
VB 2 - 1 

B - cosh (c+ {z-Vj- A) 
VB 2 - 1 

whereas 

2c+ sinh (c+ (z - L) + A) 



B2 ^ cosh (c + z + A2) 
B 2 — cosh (c + z — A 2 ) 



(42) 
(43) 



2c + sinh (c+z + A2) 
i? - cosh (c+ (z - L) + A) _B 2 - cosh (c + z + A 2 ) 
2c+ sinh (c + (z — L) — A) 2c + sinh (c + z — A2) 



i? — cosh (e+ (z — L) — A) B2 — cosh (c + z — A 2 ) 



-Ky/jmj. (44) 
-Ky/jm2 (45) 



follow from eqs. (|36|) . 

Eqs. (|42II45|) constitute a set of four equations for the four unknowns 
A,B,A 2 , and B 2 . Although their general solution is quite complicated, the 
ansatz 

B 2 =B (46) 
immediately solves (I42I43[1 as long as 

c + (z - L) + A = ± (c + z + A 2 ) and c + (z — L) — A — ± (c + z + A 2 ) (47) 

From eqs. 1)44145(1 . we find that 

A 2 = -A and z = L/2 (48) 

or else mi = m 2 = 0. After some algebra we obtain 



B 



2c+ J m2 — mi 



^V7 



cosh 



mim2 



c_i_L 



sinh a 



+ 



^±^± sinh f £±^ + ^ cosh ( C + L \ 



ra\m 2 



\ 2 



2c+ 



V 2 



cosh a 



A = tanh 1 



m2 ~ mi sinh ( ^ 

ra 17712 \ 2 



& cosh sinh 



(49) 
(50) 



as the solution to I(44l45jl . 

This yields 



* = -2 In [/? - cosh (c+/(x) + a)] - 2 In 



^n 



(51) 
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for the dilaton, where 

{x - L x > L/2 

-x L/2> x > -L/2 (52) 

x + L - L/2 > x 

is the continuous sawtooth function, with kinks at the locations of the particles. 
Substitution of JSlJ into (J23) yields 



/3-cosh(c + /(x)+a) l ° dJ 



x-z a (x°))\ (54) 



as the solution for II . 

Consider next eq. i|16|l . which can be written as 

Ng = r 7 + N |4 + K /2^2 V7 ™ a 5{ 

Proceeding as before, the general solution is 

N 0I = + Ai cosh (c+ (x + (1-2) L)+d) (55) 
c + 

The smoothness conditions at the join are the same as in eqs. (|34|l , with 'J 
replaced by -/V . These imply Ai = A3 and C\ — C3 = C. The same is true 
for conditions (|35|l . which then yield Ai = A2 and C2 = —C. The junction 
conditions (|36[) . however, are replaced by the conditions 

NM - N M = ^mxN m {z) N' Q2 {-z) - N' m {-z) = ^m 2 N m (-z) 

(56) 

where z = L/2. After addition and subtraction we find 



Ar = -gr C = a (57) 



iV = (cosh (c + /(» + a) - #) (58) 



so that 

T7 

where f(x) is the sawtooth function given in 1)521) 
Integration of fTTjl yields 



t 

where the Lfi are constants of integration and 



N11 = 7— /» - sinh ( c +f( x ) + + K * ( 59 ) 



x - L a; > L/2 

/ (x) = <( x L/2 > .x > -L/2 (60) 

x + L - L/2 > x 
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is the discontinuous sawtooth function, with discontinuities at the locations 
x = ±w of the particles. However the function Ni (x, t) must be continuous at 
these locations and smooth at x = ±L, which implies that 



fic+C+L — JTT 



-Kn 
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7 TT 

2c%(3 



smn(— h a) + smn(— a) 

smh(— h a) — smh( — a) 



0(61a) 
(61b) 



where the latter condition holds since iVi (x, t) must vanish at the locations z a 
of the particles from eq. (|21J) . 

Employing the dchnition of (3 from <|49[) we find from 1)6 la|) 

m-2 — mi . d f c+L 
sinha— sinn 

<ZT 



^Icoshf^+^Binhf^ 

UTO 2 V 2 / 2c + V 2 

which from (|5(J|) integrates to 



- cosh a— 



m 2 
mi 



= (62) 



cosh 



c+L\ nMy/ 7 ) 



2c 4 



sinh 



c+L 



~] 2 



A 2 sinli 



c+L 



= e 



where 



M 



7774TO2 



A = 



m-2 + mi 
and £ is a constant of integration. 



7712 — mi 
m,2 + mi 



(63) 
(64) 



Inverting this equation to find y/j as a function of r, we find 



Lyjr 2 - A/2 
and subsequently 



=arctanh 



2 i\i 2 



(e 2 -l)(l-AVe 2 )(r 2 -A/2)~^ 



-£ 2 (l-A 2 /£ 2 )(r 2 -A/2) 



Vr 2 - A/2 



(65) 



= 



arctanh 
2 



A- 



kM ^2 



+ (£ 2 - 1) (1 - A 2 /C 2 ) (r 2 - A/2) 



(a 2 - e) {r 2 - a/2) + 



(A-l) 2 (r 2 - A/2)- 



? 2 \ - 1 



(A + l) 2 (r 2 -A/2)- 



kMA (r - A/2) (A 2 + 1 - 2£ 2 ) (r 2 - A/2) - 



-£ A' + 1) (t' - A/2) - 



> 2 - 1) (r 2 - A/2) 



+ (e 1) (1 - A 2 /C 2 ) (r 2 - A/2) 
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from eqs. (|49I5U|) . 



cq. (fTHfl can be simplified to 



pn(«n)]- 



N[ N 



i[ln( K n)]' = 



(66) 



since iVoexp\I>/2 is independent of x. Insertion of and into the left- 
hand side (|6T)|) yields 



1 



cosh (c + x — a) — (3 



(3(3 



2 

7 TT 



((3 cosh (c + x — a) — 1) 



sinii ( <■ ... .r - n ) ( d — — 2 cosh ^— t— ^ sinh (a) 



(67) 







in the region |x| < L/2. 

The preceding equation l|67|) has no solutions unless A = 0. The problem is 
that each term has to vanish separately, ie. 



9 

7 TT 



3 2 



1 



= and 



2 

7 TT 



cosh 



sinh (a) = 



and this is not allowed from the definitions of a and (3. This peculiar situation 
suggests (within the context of our ansatz) that only the equal mass case is 
allowed as a solution to the field equations due to the consistent time-evolution 
of II. It is remarkable that the unequal mass case can be exactly integrated, 
yet the II equation cannot be solved. When this calculation is repeated for the 
c 2 = (?_ < case the same result occurs: the field equations can only be solved 
when the particles have equal mass. 



4 Two-Particle Case 

Because the total charge of the system must be zero there are no solutions for a 
single charged particle. We consider in this section the two-particle case, which 
is the simplest non-trivial charged circular self-gravitating system. 

Since there is freedom in choosing the origin we choose the origin to be 
halfway between the two particles, which are thus are located at z and —z. 
Since the total charge is zero we can assume (without loss of generality) that 
the particle at — z has charge q and the particle at z has charge —q. Therefore 
the electric field is given by the equation E — ^q(e(x + z) — e(x — z)) +E in 
this case (and consequently E 2 = (^q 2 + qE n )(e(x + z) — e(x — z)) + Eq). 

We next consider solving the field equations (|23II - H3UI) . The solution to 
equation (|2"3|) yields the result 

n = n (t)exp(*/2), (68) 
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where Ho(i) is an arbitrary function. It is useful at this stage to introduce some 
new variables 



A OU i = A — kEq, (69) 
K in = A — n(q 2 + 2E Q q + E 2 ) (70) 

so that in equations H24|) and H25[) the term kE 2 is essentially just an adjustment 
to the value of the cosmological constant A. This is similar to the result of (1+1) 
lineal gravity with charge in which the electric field becomes sort of a "local" 
cosmological constant. By introducing A out and Aj„ it is possible to write 

7 (r 2 - A/2 + nE 2 /2) = ( ^ - ^i^} "i "I" 11 h > " ) (71) 
n ' I 7( T ~ Aj„/2) = d 2 when \x\ < z J v ; 

where x is periodic with period 2L. 

We shall assume that c 2 = c\ > 0, d 2 — d 2 ^ > 0. The analysis for c 2 — c 2 _ < 
0,d 2 = d 2 _ < is fully analogous and we shall only discuss the results for this 
alternate case. 



4.1 Solving for 

In order to solve for Nq we will first integrate (|25|l across the particles to get 
the following two matching conditions across the particles for Nq. 

N' {z + e)-N' Q {z-e)='^^N {z), (72) 
N -(-z + e)-N -(-z-e) = ^lN (-z) (73) 



Then we also assume that Nq should be a continuous function, which gives 
the following matching conditions 

No(z-e) = N {z + e), (74) 
N (-z-e)=N (-z + e) (75) 

where the limit e — > is taken throughout. 

Solving Ij25(l individually in the regions between the particles is straight- 
forward and leads to a general solution for 



!Acosh(d+x) + E sinh(d + x) — \x\ < z 
+ 
Csmh(c + (L-\x\))+Dcosh(c+(L-\x\)) - ZJUaH > z 
c + 

, where A(t),E(t), C(t) and D(t) are x-independent. We can then use the four 
matching conditions across the particles and the smoothness conditions to solve 
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for the four unknowns A, E, C, D:. 

E = C = (77) 
= cosh(ri + z) r 7 ( 1 1 \ 

cosh(c+(z-L)) cosh(c + (z-L))V4 4/ 



! ^)(^ +c+d + tanh(c+(z - L)) (^-^)s ,7!)) 



+ 

where 

B = — — — cosh(d + z) + d + sinh(d + z) — c + cosh(<i + z) tanh(c + (z — L)) (80) 

or we can write D in the form 

= nm^j ( JH_\ cosh(d + z) 

2 \Bd 2 + ) cosh(c+(z-L)) 1 j 

+ tanh(c + (z - L)) coah(M) + l) _ £)) (-1 - -1) 

Note the resulting solution also has the following the symmetry properties 

N (-z + x) = N (z-x), N {~z~x) = N a (z + x) 
N' Q {-z + x) = -N' Q [z-x\ N' Q {-z-x)=-N' Q {z + x) 



4.2 Solving for \l> 

To obtain the solution for ^ we make use of an ansatz previously employed in 
the neutral TV-body case £Q 

* = -2H-N (x)^) 2 In (-^=) (82) 

where 

»-{EB<:} < 83 » 

Htli 1 ^} < 84 » 

Note that [3 ou t and /3j„ are constants in terms of x. Using this ansatz along 
with equations (|25|l and (|68|) we find 

*< = -2*° 

•* - -»#W£V~£W#V-»'-E^"C«-«o 



7V V^o/ V^c 

(nW 2 - 1) 

TV 2 /? 2 /! 2 



(«n) 2 = (Kn ) 2 exp(*, - 
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Inserting these results into 1 (24(1 leads to the equation 

(N'f - 2f 7 7V - h*N 2 - ( ' 7) y^ 1} = 

and substituting the equation for Nq into each region individually allows the 
calculation of (3 ou t and /3j„. For example when |x| > z 

(ri? 



c +Pout 
7"7 

Del 



^<(sinh^( C+ (L - \x\)) - cosh^( C+ (L - |x|))) = - ;J ;J 

or alternatively /3 ou t 
A similar calculation for \x\ < z yields 



T7 



To summarize, the solution for I24|) is given by the equation 



f -21n(- W )-21n(^|^),N 



> z 



-2 In(-JV (x)) - 2 In ( j^f^ ) , N < z 
where (3 out and /3 ira are respectively given by eqs. and 



(85) 



(86) 



(87) 



4.3 Solving for iVi 

Integrating (fTTjl yields the general solution for N\ 



(i 



2 _l_ 7 TT 

2 ~ "ST" 



7 _|_ 7'Tf j ^ _ - rl) 



x — sinh(c+(:E + L)) + _ftTi, x < —z 



sinh(d + a;) + —z < x < z 
sinh(c + (x — L)) + K%, x > z 



it A 



(88) 



To find the unknowns K\ , K 2 , and K3 in eq. 188(1 , it is convenient to use the 
constraint equations I3U|) for N\ . The constraint TV^Za) = implies that 

Ni(-z-e) = 0, N 1 (-z + e) = 
Ni(z-e) = 0, iVi(z + e)=0 

Together Ni(-z + e) = and iVi(z - e) = imply that K 2 = and 



7 ( — sinh(d + z) ) = 



(89) 
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where we have employed the following identities 



c+ 
7 — 

c+ 



7 7 TT 
2 + 



7 7 TT 
2 



(90) 



that follow from eq. (|T 1 11 . 

Equation (|89|l is a new relation between the functions d+,7 and A. Using 
N x (-z - e) = and 2Vi(,z + e) = along with the condition Ni(-L) = Ni(L) 
yields 



-K 3 = j—L 

c+ 



7 (^-i) 



t£> 



sinh(c + (z — L)) 







7^±(x + L) - ^ sinh(c+(x + L)),x < -z 



G+ - ' c + 

7 £(z - L) - ^ sinh(c + (z - L)), x > z 



7^-x — 3j— sinh(d+x), —z < x < z 



with (|30f) now being represented by the constraints 



(91) 
(92) 



where (|92f) is another new relation between the various functions. So the 
solution for Ni is 



(93) 



d + tA ■ \,fj \ n 
— — z — sinn(a + z) = I) 

d + d + 

— (z-L)-—smh(c+(z-L)) = 
c+ c+ 



(94) 
(95) 



4.4 Constraint equations 

We must now verify that the constraint equations Ij29(l . (|94(I . H95(I can all be 
satisfied. It is useful to introduce a new parameter 



n = = y(g + 2£ ) 



so that = c 2 + + 771, or 



1 

4 



7n 



(96) 
(97) 



The constraint (|29|l can then be written as 







<i+ tanh(d + z 



2^7 



. ft/ft . 9 



77Z 

2 — d + tanh/d+z) tanh(c + (z — L)) 



+c+ tanh(c+(z — L)) 



27271 / Km 



rnncZ 



((=) 



V7 ^tm s 2 
me V 2 j 



(98) 
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where the calculational details are given in the appendix. Rather than first 
attempting to solve this equation in full generality, we instead consider specific 
values of n. 

The simplest case to examine is when n — 0, or q = —2Eq. eq. I|98|) reduces 

to 

= c +^~^ — (tanh(c + z) + tanh(c + (z — L))) 

which in turn implies that z = ^, and so the particles must be located at -j 
and — The electric field is therefore given by E — hq(e(x + z) — e(x — z)— 
1) or more explicitly 

E = -iq, |.t| > z^q, \x\ < z 

Note that E 2 — which is constant, so the magnitude of the field doesn't 
change with position. 

We next rewrite l|94j) an d l|95Jl in a more convenient form. After some 
manipulation, 194(1 can be written as 

d 2 + d + z ( — ^- — cosh(d + z) + d + sinh(<i + z) — c + cosh(d+z) tanh(c + (z — L)) 

JT r ^ Km V^ + p tanh(c+(z - L))^j sinh(d+z) = (99) 

and (|95|l can be written as 



c + (z — L) — t ( ) ( — + — tanh(c + (z — L)) ) cosh(d + z) tanh(c + (z — L)) 



2 2 I tanh(c + (z — L)) = 



(100) 



Since n — 0, c + = d + . Recalling that z = k, we have 



c+c+- ^ — - — cosh(c+-) + 2c+sinh(c + -) J - jtt — - — sinh(c+-) = 

(101) 

for 199|) . It is straightforward to show that (|100l) reduces to the same equation, 
cq. I|101|) can be integrated to yield 

-L. Km^f 



cosh(V7 {r 2 - kout/2)-) + sinh(V7 (r 2 - A out /2)-) = £ 

z 4c+ 2 

where £ is a constant of integration. It is then possible to solve for 



'7= arctan fa 5 — 5 v( T A out /2) 

- A onf /2 \ ^f +e{T 2_ Aout/2) VI out/ , 

(102) 
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One might expect that a solution of for general n would provide infor- 
mation on the locations of the particles. However due to the extra terms present 
when n ^ 0, it does not appear that either of these constraint equations is a 
total derivative. For example, consider setting n = (-^) 2 so as to simplify the 
equation as much as possible. In this case eq. <|98H reduces to 

= d + tanh(rf + z)Km^/7 1 + tanh(c+(z — L)) 

This implies either 

d + z = 

or 

tanh(c+(z — L)) — 



We can only satisfy the first equation by setting z = 0. However this corresponds 
to both particles at the origin, i.e. a single particle at the origin with mass 2m 
and no charge; this was previously solved |17| . 

Turning to the second equation, since z is independent of time we can regard 
this equation as yielding a solution for 7 as a function of r. A problem occurs 
because there are more constraint equations that need to be satisfied, and in the 
uncharged case the other constraint equations are the equations that determine 
the time dependence of 7. 

In this case the other constraint equations cannot be satisfied. This can be 
most easily seen by applying this result to l|100|) . Substituting tanh(c+(z— L)) = 
~~ Km^/f an< ^ n = (' i T i ) 2 m t° CUnt yields the following result 



c +(*- L ) + ^(^T-) = ° 



Now d\ = c\ + 7 n = c\ + 7 (^f) 2 = 4(1 + (^5) 2 ) 



so 



. = 1 nm^j 7 rf 1 

C+ z - L 2c+ 4 1 + ( "W7 )2 1 ' 

But tanh(c + (z — L)) = — K ^/7 gives us a solution for c+ i.e. 



^tanh- (^±=) = ^tanh- 



z — L \ Krriy/j J z — L 

which implies that 



c + = -- L- i ^V77I! (1Q4) 
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Comparing eqs. I|1U3|) and (|104f) implies that KT ^~^ = , ie. c + = which 

contradicts our initial assumption that c+ > 0. The same problem occurs 
when the case c 2 = cL < when looking at the case n — ( ! ^) 2 ^ 

Another simplification occurs if n = — (^f 1 ) 2 , Similar manipulations yield 

tanhfa+zl = 

which is of the same form as the previous result. Again this is a problem 
because the other constraint equations are not satisfied. While these results 
don't rule out solutions for all n =/= , they illustrate the problems in finding 
them. In fact intuitively one might think that n — — (^) 2 and n — (^f 1 ) 2 
provide bounds on n. 

We have been unsuccessful in finding any further ansatze that could solve 
the constraint equations. Further progress might be attained by attempting 
to solve the constraint equations numerically for different values of n to see if 
solutions exist for n^O. However our initial attempts at solving the equations 
numerically have uncovered further difficulties: the problem is that the system 
appears to be overconstrained when n ^ 0. When n = then and (|100|) 
both reduce to the same constraint, but this not the case for arbitrary n. 

We turn now to consideration of the n = case. 



4.5 The n = solution 

When n = the previous solutions of Nq,^,Ni simplify and the constraint 
equations derived from i|29|) and (|30[l are satisfied. The solutions can be written 
as follows 



II = n (t)exp(tf/2), (105) 

T7 T7 

T 

+ 



N o = 7T^cosh(c + /(x))-^ (106) 
pc_|_ 



where 



and 



^ =cosh(c 4 )+ ^ sinh(c 4 ) 



f(x) = {\ x \, \x\ <^L- \x\, \x\ > ^} (107) 

is the saw-tooth function that has magnitude zero halfway between the particles 
and magnitude k at the locations of the particles. We also have 

+ L) - sinh(c + (z + L)),x<-% 
•V, = { J^x-^siuh(c + x),-±<x<% } (His) 
!%{x - L) - smh(c + (x - L)), x>\ 
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Notice that eq. (|105fl an d (|109fl together imply 



(3- cosh(c + f(x)) 

and solving (|27|l we obtain 



«n=/ c +\f?„ (no) 



Using the solution (|102|l or the metric 

2 ^v/^ + (e 2 -l)(r 2 -A e /2)-^ 

\/7 = : arctanh — =— = V (r 2 - A e /2) 

V/ L-y/ t 2 - A e /2 I ^+^ (r 2_ Ae/2) VV 



(111) 



allows us to rewrite 



P=±^^- + { e-l)(T>-A e /2) (112) 

These equations are the solution for when c 2 = c+ > but solutions also 
exist when c 2 = c 2 _ < which can be obtained by substituting c + — * ic_ into 
the above results. The solution for the metric function then becomes 



V7 = 



H HI 

4 



arctan — 5 — 5 J (A e /2 — 

\ ^ - e (A e /2 - r 2 ) 

(H3) 

where fc is an integer and |£| < \ fl + as a consequence of the periodicity 

of Ni(x, t). No solutions exist for c 2 = 0. 

Note that in either case the two particles remain in equilibrium in the sense 
that the ratio of the proper lengths between them as computed from either 
the rightward or the leftward side of the particle at z = is constant. These 
solutions are very similar to those obtained in the neutral system the key 
difference being that cosmological constant has been replaced by A e = A — kE 2 
. The electric field is constant in magnitude but alternates in sign between each 
the particles. This result can be easily generalized to N particles. 

5 Charged N-Body Solution 

The only difference between eqs. (|23fl - l|30|) and their counterparts in the neutral 
self-gravitating system is the presence of E 2 in equations l)24 [l .l|25 [l and the 
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term J dxNoy/jEjj^- in equation 1|29[1 . Using the ansatz E 2 — constant, it is 
possible to eliminate E 2 in equations l|24 |l .l|25 |l by introducing A e = A — kE 2 . 
Furthermore, we have 

f , dE , dE , ?MZ a +e f Za+e , dE 

dxE = / dxE — =(E 2 )\ a+ - dxE — 

J dz a J Za _ e dx ,z « e J Za _ e dx 

[ Za+e dxE ^ = i(s 2 )r» +6 = o, 

Jz a -e dx 2 1 

when E 2 is constant. This implies that J dxN 0v /jE-^M- — 0, since J dxN 0y /^E^j- = 

So an equilibrium solution for N particles with charge distributed evenly 
around the circle exists when E 2 is constant. This yields for c\ > 



\Fi = ; = =arctanh f WN 2 , f2 — J t 2 - A, /'2 

Wr 2 -A e /2 \ ^l +e{T 2_ Ae/2) 



(114) 



where M = mN is the total mass of the system. We also have 



* = -2 In 09 - coahfo /(*))) - 2 ln( ( * } ) (115) 

c+Vf 3 - 1 

No = ^(cosh(c + /(a:))-/3) (116) 
c_j_p 

2 • 

Ni = 7— h(x) - X^- smh(c + h(x)) (117) 

where f(x) is the saw-tooth function that peaks with a value of L/N (i.e. 
f(z a ) = L/N) at the particle locations and has zero magnitude half-way be- 
tween the particles (i.e. f{ Za+ * a+1 ) = so /( * tt+ * a+1 + y) = |y| for < L/N). 
The function /i(x) is a function that has a jump discontinuity across the par- 
ticles, and between two particles linearly changes from —L/N to L/N (i.e. 
h( Za+ * a+1 +y) = y, for \y\ < L/N). Notice that \h(x)\ = f(x) and h'{x) = 1 
between the particles. 

The functions (3 and Uq are given by 



AN k 2 M 2 



, , K 2 M 2 



no = ^^iy^^ + ^-A^J (119) 
where £ and £ are integration constants. Then 



cosh(c+/(x)) - (3 
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It is possible to write the functions for H!,Nq,Ni in a different form by using 
the following identity 



cosh(c+/(x)) = , N ' V cosh(c + (|x - z a \ - L)) (121) 
sinnic+L) ' 

x ' a— 1 

whose proof is given in the appendix. Then for finding N\ it is possible to 
integrate l|26[l to find the following expression 

c ^ 2 tt sinh( c+ ^) ^ 
Ni = -y—x g— - N ^ e ( x ~ z a) [sinh(c+(|a; - z a \ - L)) + sinh(c + L)] 

As with the two particle case there are solutions for c 2 — c 2 _ < 0. These are 
obtained upon replacing c + — * ic_ in the above results. The solution for the 
metric function 7 now becomes 



r N 



^ + (l-a(A e /2-T 2 )-if 



arctan — 5—= \/A e /2 — r 2 + k-n 

1 ^_ C 2(A e /2-r2) V 1 

(122) 



where k is an integer and the integration constant |£| < y 1 + sa^Iv 2 as a con- 
sequence of the periodicity of N\(x, t). 

So an equilibrium solution exists for N charged bodies when when E 2 is 
constant. The charges of the particles alternate in sign around the circle; since 
the total charge vanishes we have a solution only when N is an even number. 

If the total number of particles is odd then some particles must remain 
uncharged. The simplest example is the 3-particle case, one particle has charge 
q, one has charge —q, and one particle has zero charge then it is possible to 
choose Eq so that E 2 = const. 

The previous results can be approximated for the case when we take the 
limit N — > 00 but we keep the total mass of the ring M fixed . 

For the case that c\ > then the solutions can be simplified to: 



x /7 = ; = =arctanh ( - ^ — 

V ' ZVt 2 - A e /2 I t 



c 2 , = — -arctanh 



L 2 \ £ 

AN 



= —^(e-l)(r 2 -A e /2) 



n, 



V^C 2 (r±v/(r 2 -A e /2)) 
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* = -2 In 09 - cosh(c+/(aO)) - 2 ln( fd ^ff_ ) (123) 



iV = (cosh(c + /(x)) - 0) (124) 



2 • 

JVi = j^h{x) - X^- siuh{c+h(x)) (125) 

where f(x) is the saw-tooth function that peaks with a value of L/N (i.e. 
f(z a ) — L/N) at the particle locations and has zero magnitude half-way be- 
tween the particles (i.e. /( z "~^ a+1 ) = ). 

For the case <?_ < the solution for the metric function 7 becomes 



N 



VAe/2 - T 2 



arctan 



fc7T 



c_ = 



L 2 



arctan 



fc7T 



where fe is an integer and the integration constant |£| < 1 



6 Summary of N-body Solutions 

The 2-body equilibrium solution we found is simply a special case of the N-body 
equilibrium solution that occurs when E 2 ^constant. For convenience here is 
a table summarizing the results for c 2 + > 0: 



variable 


solution 


n 


n (t)cx P (*/2) 




-21n(/3- cosh(c + /(x))) - 21n(-^_) 




^ (cosh(c + /(x)) - /?) 




l C 4Hx)-^ S mh(c + h(x)) 


n 


V7C 2 (r ± v/^T-D + ^ 


a 


f^ + K 2 -l)(^-V2) 


4 


7(r 2 -A e /2) 




, N arctanhf W^+f"^-^-** ^^Y^ 


L^r*-h e /2 y ^ +? 2 (T 2_ Ae/2) V e/ J 


m 


f(^ ±1 + y) = \y\ for y| < L/N 


h(x) 


h(^±l+y)=y,for \y\ < L/N 


A e 


A-kE 2 
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For c 2 = c 2 _ < replace c + — > ic- and the metric is now 



variable 


solution 


N 


--^(cos( C _/(z))-/J) 


Ni 


7^) + ^sin(c_^)) 




w 


arctan ( ^^^/^"w JEJTl*^ 


+ A:7r 






y «^-^(A«/a-T») V e/ y 




iei<>/i+eft 



where fc is an integer. 



7 Discussion 

We have obtained an exact equilibrium solution for N charged particles of equal 
mass in circular gravity. Our notion of equilibrium has been taken to be that 
of a static equilibrium: namely that the ratio of the proper lengths between any 
two pairs of particles is constant. On any t ^constant hypersurface this implies 
that the values of the z a are constant. 

Our solution is valid only when the electric field is of constant magnitude. 
The electric field will be constant when the number of charged particles is even 
and these particles must have the same magnitude charge alternating in sign 
between positive and negative around the circle. Although restrictive, there 
is still a lot of freedom in choosing the arrangement of the charged particles 
situated amongst any number of neutral particles that are present between them. 
For example it is possible to have a system of N particles with only two charged 
- then there would be N — 1 ways to arrange the charged particles within the 
system. 

Apart from these distinctions, our solution is otherwise the same as that 
obtained in the electrically neutral case pp. Our solutions almost certainly 
describe an unstable equilibrium, since the masses are all equal and the particles 
are evenly separated. In fact an A-body stable equilibrium does exist within 
the context of our restrictions, but corresponds to all particles on top of each 
other, which is simply the one body solution. 

We were unable to obtain the most general analytic solutions to the full set 
of equations. An inspection of the constraint equations suggests that they can 
only be simultaneously satisfied under very stringent circumstances. Indeed, 
the situation with two bodies of unequal mass is quite perplexing. One might 
expect that a solution of unstable equilibrium would exist for two bodies located 
equidistant from one another on the circle. However our rather general ansatz 
(|46|) yields a consistent solution to all equations except for that which governs 
the time-evolution of II (eq. UlSfl ): we can find solutions to this last equation 
only when the masses are equal. Similar stringent constraints hold in the 
charged cases, where again we were able to find solutions only for equal masses. 
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It would appear that the most promising way to make further progress would 
be to solve the equations of (1 + 1) circular gravity numerically. Lineal self- 
gravitating systems have been studied numerically and interesting results have 
been obtained in both the 2-body JS] and 3-body [2] cases. Currently there 
have been no numerical models of (1 + 1) circular gravity. It would be interesting 
to numerically model two particles on a circle and to see how the space-time itself 
behaves differently when compared to the behavior of two-body equilibrium 
solution. It would also be interesting to numerically model the system of 
particles when the masses are unequal or when charge has been introduced or 
both. 
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8 Appendix 

8.1 Derivation of Equation (98) 

In examining equation i|29|) notice that N' (z a ) needs to be defined since N' (x) 
is discontinuous across the particles. This has been seen before in equation l|73|l 
So N' (z a ) is defined as 

,,,, s_ l; N' (z a -e) + N' (z a + e) 

N {z a ) = km (126) 

e— »0 Z 

Now examining the / dxNo^E^j- term of equation it is useful to rewrite this 
term as 

J dxN 0v /jE^- = N 0v /j J dx ^ e b e{x - z b ) + E^j e a 5(x - z a ) 

= No(z a )^/e a ^ ^ e b e(z a - z b ) + Eg j 
Returning to the two particle example, this implies that at x = —z 



dxN 0y /jE 



dE 

d(-z) 



N (~z)^-q(q + 2E ) 

N (-z)Sr- 

K 



and at x = +z 



I 



dxN y^/E 



dE 

dz 



-N (z)^y-q(q + 2E ) 
-N (z)^y- 

K 
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Note that because of the symmetry properties of No we have Nq(—z) = 
-N' (z) and / dxN s/jE-J^ = - j dxN Qy ^fE^. So JgjJ is the same equa- 
tion at z and —z and therefore it is only a single constraint equation. Explicitly 
Nq(z) is given by 

1 ( i~y n *| 

Nq(z) — — < A [d+ sinh(c?+z) + c+ cosh(e?+z) tanh(c+(.z — L))] H tan(c+(z — L)) > 

2 [ c +"+ J 



where 



and so 



7-7 _ AS 



. . Acosh((i + z) 1 — , . , , 

^oC*) = , = LJ - J r — ^d+tanh(rf + z) 



2(^ + ^tanh( C+ (.-i))) 

^ — c+ tanh(c+(z — ij) 1 

771 



(127) 

C4 





1 


2 












> 



The constraint equation l|29[l implies that 



KM = 



niH 



T7 

^4cosh(d + z) — -p- 



or alternatively 



Acos\\(d+z) 



( 2771 

2 (^^ + ^tanh(c + (z-i))) \V^ mK 
2jn 



c + tanh(c + (z — L)) 



771 
75" 



d+ tanh(d_|_z) 



(128) 



Since the constraint (|128|l must be equivalent to the definition of N' {z) in 
eq. I|127|l we have 



= d + tai\h.{dj r z) 



2^7 



. Kill . n 



77X 

2 — d + tanh(<i + ;z) tanh(c + (z — L)) 



+c + tanh(c + (z — L)) 



2yzn /.nm 2 



((f) 2 



- n + 



2^/7 / Km 



iriK 



((f) 2 -) 



129) 



which is the constraint equation that results from (|29|) . 
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8.2 Proof of Identity (121) 

Now returning to the TV-Body equilibrium solution when E 2 = const an inter- 
esting identity that was mentioned without proof was eq. I|121|) . This identity 
is useful because it allows us to write the functions ^f,N ,Ni in two different 
forms: either as a single function or as a sum of N functions. Without this 
identity it would not be apparent that both forms of the equations are in fact 
equivalent. 

We now verify eq. Q121fl , We begin by noting that both cosh(c+/(x)) and 
sinh(c+£) S^=i cosh(c + (\x — z a | — L)) arc periodic between the particles. This is 
equivalent to saying that both functions have a period of For cosh(c+/(x)) 
this can be seen because the sawtooth function f(x) peaks at the particles, 
yielding f(x+^k) = f(x). In fact there are two values for each \x — z a \ because 
this system is located on a circle. Thus there are two possible directions in 
which one can measure the distance between x and z a . If \x — z a \ is the distance 
between them one way around the circle then 2L — \x — z a \ is the distance as 
measured via the other way. Note that regardless of how \x — z a \ is measured 
the functions cosh(c+(|a; — z a \ — L)) are well defined. This is because 

cosh(c + (|ir — z a \ — L)) = cosh(c + ([2L — \x — z a \] — L)) 
= cosh(c + (L - \x - z a \)) 

Since the particles z a are evenly spaced with spacing ^ it follows that 
^2a=i cosh(c + (|x — z a \ — L)) must be periodic with period jjf. This is because 
a shift by ^ can be accounted for by relabeling the particles i.e. 

N 2L N 2L 

^cosh(c + (|x+ — - z a \ - L)) = ^cosh(c + (|a;- (z a - — )| - L)) 

a—l a—1 

JV 



= cosh(c + (|a: - z b \ - L)) 



6=1 



Therefore because both cosh(c + /(a;)) and Xh(c+L) J2a=i cosh(c + (|a; — z a \ — 
L)) are periodic between the particles it is only necessary to show that these 
functions are equivalent between two particles in order to prove that they are 
equivalent on the entire circle. To prove the equivalence it is necessary to 
separate the cases when N is even and when N is odd and prove these cases 
independently. 

1. Even N 

For the case when N is even the origin is located halfway between two 
particles, and so the particles are located at J L for k = — + 1, — y + 2, 

y . We will show equivalence of cosh(c+/(ir)) and ^^"^ S a =i cosh(c + (|a; — 



2G 



z a \ — L)) when ^ < x < j^. Recall that the sawtooth function f(x) is defined 
to have value at the particles (i.e. /( 2A ^ 1 £) = j^) and the value halfway 
between the particles (i.e. f(^-L) = 0). So this means 

fix) = \x\ for — \- < x < ^~ 
j v J ii N - - N 

Now 

cosh(c + |x|) = cosh(c + x) 

.-. cosh(c + ,f (x)) = cosh(c + a;) for < x < 

and so it is only required to show that J2^=i cosh(c + (|a;— z a \— L)) = sinh (^+^) cosh(c + x) 

sinh(-fc-) 

This is straightforward 



N 



^2 cosh(c+(|a; - z a \ - L)) 



a=l 



fe=i ' fc=-f+i 

f o 

= J2 cos Hc+(^—L-x-L))+ cosh(c+(a; —L-L)) 

fe=i " fc=-f+i 

* gfc-CJV+l) , 2fc-(iV + l) r , 
= 2^cosh(c + ( i-a:)) +cosh(c+( L + xj) 

fe=i 



cosh(c+(|a; - z a \ - L)) = ^ 2cosh(c + a;) cosh(^(2fc — (TV + 1))) (130) 

o=l fe=l 
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and note 



^ cosh (£±^ (2 fc_(7V + l))) 



k=l 



fe=l 



!+^(2fc-(JV+l)) + e -^(2fe-(iV+l)) 



E 

k=l 



sinh(^(2fc - N)) - sinh(^(2fc - (JV + 2))) 



e n — e jv 



1 



2sinh(^) 

sinh(c + L) 
2sinh(^) 



2 sinh(^(2fc - N)) - £ sinh(^(2(fc - 1) - AT)) 



fc=i 



JV 



Jc=l 



So we have 



E, ,c + L . , sinh(c+L) 
cosh(-±-(2fc-(JV + l)))- 



fc=i 



2sinh(^) 



and by substituting (|131|) into 1|13U|) we obtain 



(131) 



JV 



V cosh(c+(|x - z a \ - L)) = bm ^^ cosh(c + a;) 
Hence for an even number of particles 



cosh(c+/(x)) = \ w W n Vcosh(c + (|x-z a | - L)) 
sinn(c+L) ' 

v ' a— 1 



2. N is odd 



When the number of particles is odd the origin is located to be at the location 
of one of the particles, so the particles are located at for k = —^r^, —^^- + 
1, ^y^- This time we will show the equivalence of cosh(c + /(a;)) and 

71\Tl) cosh (MN - z *\ - L )) in the re S ion < x < f . 

So 

/(0) = /(f) = | 



L x JV iul JV - x - JV 



2N 



Note 



cos h( c +(^ - x)) = cosh(c + (.x - — )) 

L 2L 
cosh(/(x)) = cosh(c + (.x — — )) for < x < — 



By a similar calculation to the even case we can show that 



a=l 



cosh(c + (|a; — 



z a \-L)) = ^^-cosh(c + (x-§)). 



Hence when N is odd 



cosh(c+/(x)) 




a=l 



^cosh(c+(| 



N 



x - z a \ - L)) 
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